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1  INTRODUCTION 

The  Air  Force  Geophysics  Laboratory's  LOWTRAN  Model  (Kneizys  et  al., 
1980#  1983)  is  used  to  calculate  atmospheric  extinction  in  the  United  States 
Air  Force  Tactical  Decision  Aids.  LOWTRAN  uses  a  set  of  model  aerosol  pro¬ 
files  to  represent  extinction  properties  of  various  types  of  atmospheric 
aerosols.  A  brief  description  of  these  models  follows.  A  more  complete 
description  is  given  by  Shettle  and  Fenn  (1979) .  LOWTRAN  aerosol  models  in¬ 
clude  the  following: 

-  FOC  1  characterizes  very  wet  aerosols  found  in  very  dense  fogs, 

-  FOG  2  describes  aerosol  properties  in  less  dense  fogs  such  as 

radiation  fog. 

The  present  report  describes  other  non-fog  types  of  aerosols. that  exist 
at  relative  humidities  below  99  percent.  The  four  LOWTRAN  aerosol  models  that 
characterize  the  non-fog  aerosols  are  the  rural,  urban,  maritime  and  tropo¬ 
spheric  models. 

(1)  The  rural  model  describes  the  boundary-layer  background  aerosol 
found  in  continental  airmasses.  This  aerosol  is  composed  principally  of 
sub-micrometer  soil  particles  and  products  of  certain  natural  and  man-made 
processes  occurring  in  continental  regions*  The  particle  size  distribution 
in  the  rural  model  represents  a  concentration  of  particles  in  the  size  range 
below  1  pm  and  contains  two  modes,  one  at  about  0.03  pm  ard  a  second  mode  at 
0.5  pm.  The  number  concentration  in  the  second  mode  is  very  low  (about  1.25  x 
10  times  that  in  the  first  mode). 

(2)  The  urban  model  describes  aerosols  that  include  significant 
concentrations  of  particles  produced  in  urban  and  industrial  complexes.  The 
urban  model  size  distribution  is  identical  to  that  of  the  rural  model.  The 
chemical  makeup  of  the  aerosol  is  that  of  the  rural  aerosol  and  an  added 
component  representing  soot-like  aerosols. 

(3)  The  maritime  aerosol  model  characterizes  extinction  by  the 
continental  background  aerosol  and  by  sea-salt  particles  found  in  maritime 
airmasses.  This  aerosol  is  composed  of  both  sub-micrometer  and  micrometer 
(up  to  20  pm)  sized  particles.  The  size  distribution  contains  the  smaller  of 
the  two  modes  in  the  rural  model  and  a  second  mode  representing  the  relatively 
large  sea-salt  aerosols. 


7 


the  two  inodes  in  the  rural  model  and  a  second  mode  representing  the  relatively 
large  sea-salt  aerosols. 

(4)  The  tropospheric  model  characterizes  the  aerosols  in  very  clean 
atmospheres  and  in  the  free  atmosphere  above  the  boundary  layer.  Its  size 
distribution  contains  only  the  smaller  of  the  two  modes  in  the  rural  model; 
chemical  composition  is  that  of  the  rural  model. 

Table  1  presents  a  comparison  of  8-12  ym  aerosol  extinction  coeffi¬ 
cients  from  LOWTRAN  V  for  three  of  the  aerosol  models,  for  a  range  of  visibili¬ 
ties  and  relative  humidities.  This  table  also  presents  for  comparison  the 
molecular  extinction  coefficients  at  8-12  ym  for  the  same  relative  humidity 
values  at  two  different  temperatures.  Fig.  1  presents  the  8-12  ym  transmis¬ 
sion  on  a  2  km  path  under  the  indicated  conditions  for  the  three  models  as  a 
function  of  relative  humidity.  The  following  features  of  these  data  high¬ 
light  the  importance  of  aerosol  model  selection: 

(1)  The  aerosol  extinction  at  8-12  ym  for  the  rural  and  urban  models 
is  small  except  at  low  visibilities.  With  moderate  to  high  relative  humidity 
and  visibility  greater  than  4  km,  aerosol  extinction  in  continental  aerosols  is 
overshadowed  by  molecular  extinction. 

(2)  The  relative  humidity  response  of  the  rural  and  urban  models  is 
weak.  In  comparison  the  maritime  model  shows  a  strong  response  at  relative, 
humidities  above  70%. 

(3)  Aerosol  extinction  in  the  maritime  model  can  be  relatively 
high  even  with  good  visibilities. 

Table  2  and  Fig.  r  present  corresponding  data  frbm  three  LOWTRAN  VI  models 
at  1.06  ym.  Note  the  difference  in  path  length  and  visibility  between  Figs. 

1  and  2.  Molecular  extinction  is  negligible  at  1.06  ym  and  is  omitted  from 
this  data  set.  Note  the  following  features: 

(1)  Relative  humidity  response  of  each  of  the  models  is  low.  The 
difference  in  relative  humidity  sensitivity  is  particularly  striking  for  the 
maritime  model  in  comparison  to  that  shown  in  Fig.  1. 

(2)  The  differences  in  1.06  ym  extinction  among  the  models  are 
greater  at  all  relative  humidities  than  the  differences  at  8-12  ym. 

(3)  Aerosol  extinction  is  much  stronger  at  1.06  ym  than  at  8-12  ym. 

X  . 

\  •  ' 
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Fig.  2.  1.06  pm  transmission  as  a  function  of  relative  humidity  for  three 

LOWTRAN 

Clearly,  selection  of  the  most  appropriate  aerosol  model  plays  a 
role  in  estimating  atmospheric  transmission  in  the  infrared  spectral  region. 
The  Aerosol  Model  Selection  Algorithm  described  in  Section  7  of  this,  report 
was  developed  to  give  the  users  of  Tactical  Decision  Aids  guidelines  in  model 
selection  when  estimating  aerosol  extinction  in  the  boundary  layer,  tte  based 
the  algorithm  on  published  scientific  literature  in  a  number  of  pertinent  sub¬ 
ject  areas.  This  paper  summarizes  the  scientific  data  used  in  developing  the 
algorithm  and  presents  the  underlying  rationale.' 

Section  2  discusses  methods  that  are  commonly  used  to  characterize’ 
atmospheric  aerosols.  Section  3  presents  results  of  three  aerosol  studies 
that  describe  important  properties  of  aerosols  and  illustrate  their  spatial 
and  temporal  variability.  Section  4  discusses  the  principal  processes  that 


influence  the  spatial  and  temporal  variability  of  aerosol  properties, 
tion  5  discusses  factors  that  determine  aerosol  extinction  efficiency. 
Section  6  discusses  relative  humidity  effects.  Section  7  presents  the 
Aerosol  Mbdel  Selection  Algorithm. 
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2  CHARACTERIZATION  OF  ATMOSPHERIC  AEROSOLS 

Three  properties  are  important  in  characterizing  atmospheric  aerosols 
in  the  context  ot  radiation  extinction.  They  are  the  total  number  concentra¬ 
tion,  the  particle  size  distribution,  and  chemical  composition  of  the  constit¬ 
uent  ..articles  in  various  size  ranges. 

The  total  masher  concentration  is  simply  the  total  number  of  particles 
summed  over  the  entire  range  of  particle  sizes.  The  total  number  concentra¬ 
tion  is  a  convenient  parameter  to  represent  variations  in  aerosol  concentra¬ 
tion  when  the  size  distribution  and  chemical  composition  cure  specified* 

The  particle  size  distributions  in  atmospheric  aerosols  have  been 
found  characteristically  to  contain  three  modes.  The  three  modes  are  named 
the  transient  nuclei  mode,  the  accumulation  mode,  and  the  coarse  particle 
mode.  Particle  size  values  that  define  the  limits  of  the  range  of  mode 
diameter  for  each  mode  vary  slightly  by  author.  We  adopt  the  following 
definitions  of  Willeke  and  Whitby  (1975) . 

Diameter  Mode 

<  .06  pm  Transient  nuclei 

.06  pm  -  2.0  via  Accumulation 

>  2.0  via  Coarse  particle 

The  transient  nuclei  mode  contains  the  Aitken  nuclei,  which  are  com¬ 
posed  of  products  of  both  natural  and  man-made  processes.  Examples  are 
certain  organic  compounds  released  into  the  atmosphere  by  plants,  and 
sulphate  and  nitrate  compounds  formed  either  directly  during  fossil  fuel 
combustion  or  by  chemical  transformation  of  direct  combustion  products. 

The  nuclei  mode  includes  certain  particles  that  act  as  cloud  condensation 
nuclei.  Most  particles  in  the  nuclei  mode. are  transient,  because  they  tend 
to  form  larger  particles  by  coagulation.  The  role  of  aerosols  in  the  nuclei 
node  in  radiation  extinction  is  very  small  for  visible  radiation  and  com¬ 
pletely  negligible  for  infrared  radiation,  .heir  small  extinction  cross- 
section  (see  Section  6)  renders  them  ineffective  as  infrared  attenuators. 

The  accumulation  mode  is  so  named  for  two  reasons*  (1)  many  of  the 
constituents  in  this  mode  are  products  of  coagulation  of  pc  .-tides  in  the 
transient  nuclei  mode}  and  (2)  particles  in  this  mode  tend  to  at.  emulate  in 
a  way  that  results  in  atmospheric,  residence  times  from  several  days  to  several 
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weeks.  The  small  size  of  these  particles  implies  low  sedimentation  rates;  at 
the  same  time,  precipitation  washout  of  these  particles  is  less  efficient  than 
for  larger  particles  (see  Section  4) .  Particles  in  this  mode  are  principal 
attenuators  of  visible  and  near  infrared  radiation.  Their  role  in  middle 
and  far  infrared  extinction  is  much  less. 

The  coarse  particle  mode  often  contains  particles  as  large  as  20  um 
in  significant  number  concentrations.  This  mode  includes  fly-ash  particles 
from  fossil  fuel  combustion,  airborne  dust  and  sea-salt  particles.  Residence 
times  for  particles  in  this  mode  are  short  compared  to  those  in  the  accumula¬ 
tion  mode  because  of  high  sedimentation  rates  and  efficient  washout  by 
precipitation.  Particles  in  the  10-20  ym  range  have  residence  times  on  the 
order  of  hours,  while  the  smaller  particles  (2-10  ym)  typically  have  resi¬ 
dence  times  of  1-3  days.  Particles  in  this  mode  are  the  dominant  aerosol 
attenuators  of  middle,  and  far  infrared  radiation. 

The  distributions  of  particles  sizes  within  each  mode  typically 
follow  a  leg-normal  distribution  quite  closely  (Shettle  and  Penn,  1979) . 

The  parameters  required  to  specify  the  size  distribution  for  each  mode  with 
good  accuracy  are  the  mode  radius,  the  geometric  standard  deviation  of  par¬ 
ticle  size  about  the  mode  radius  and  the  total  number  concentration  in  the 
mode.  Specification  of  these  parameters  for  each  of  the  three  modes  enables 
a  very  accurate  definition  of  the  entire  size  distribution  for  most  aerosols. 

Authors  sometimes  describe  the  aerosol  size  distribution  by  esti¬ 
mating  aerosol  particle  surface  area  or  aerosol  particle  volume  (or  mass) 
as  a  function  of  particle  size.  The  surface  area  distribution  is  useful  in 
a  study  of  radiation  extinction  because  surface  area  makes  more  apparent  the 
role  of  particle  size  in  extinction  at  various  wavelengths. 

Pigs.  3  and  4  illustrate  the  number  distribution,  surface  area 
distribution,  and  volume  distribution  for  the  mean  aerosol  condition  during 
October  1971  at  the  Denver  City  Maintenance  Yard  (Will eke  and  Whitby,  1975). 
The  modes  in  the  distribution  are  not  obvious  in  Pig.  3.  They  are  more 
readily  seen  in  the  distribution  of  surface  area  and  volume  in  Pig.  4. 

Variations  in  the  strength  of  the  individual  modes  are  useful  for 
the  present  purpose  because  they  often  directly  reflect  the  variation  in 
aerosol  extinction  in  certain  wavelength  regions. 
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Fig.  3.  The  grand  average  of  the  number  density  distribution  from  October 
1971  measurements  of  Denver's  City  Maintenance  Yard.  D  is  the 
particle  diameter  (Willeke  and  Whitby,  1975).  p 
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Fig.  4.  Normalized  frequency  plots  of  number  density,  surface  area  and 

volume  distributions  for  the  grand  average  October  1971  measure¬ 
ments  of  Denver's  City  Maintenance  Yard.  Note  the  bimodal  dis¬ 
tribution  of  volume  and  the  fact  that  each  curve  shows  features 
of  the  distribution  not  shown  by  the  other  plots  (Willeke  and 
Whitby,  1975). 
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The  chemical  composition  of  the  aerosol  is  important  in  radiative 
transfer  processes  in  three  ways.  The  first  is  the  dependence  of  the  complex 
refractive  index  on  this  property.  Secondly  the  chemical  composition  is 
directly  related  to  the  sources  and  mechanisms  for  production  of  the  aerosol. 
These  are  important  in  understanding  the  observed  temporal  and  spatial  dis¬ 
tributions  and  the  particle  size  distributions  of  various  types  of  aerosols. 
Thirdly,  the  chemical  composition  determines  the  response  of  the  aerosol  to 
changes  in  relative  humidity. 

The  chemistry  of  aerosols  may  be  classified  in  three  ways,  namely, 
soluble  or  insoluble,  organic  or  inorganic,  and  by  source,  particularly 
whether  the  source  is  natural  or  related  to  human  activity.  Solubility  is 
perhaps  the  most  important  of  these  because  it  determines  the  aerosol's 
response  to  changes  in  relative  humidity.  Examples  of  soluble  aerosol  com¬ 
pounds  are  ammonium  sulphate,  nitric  acid  and  nitrate  salts,  sulfuric  acid 
and  sodium  chloride.  Insoluble  aerosol  compounds  include  elemental  carbon 
(soot) ,  soil-derived  aerosols  such  as  sub-micrometer  clay  particles  and  much 
larger  particles  due  to  wind-blown  dust,  and  other  aerosols  due  to  human 
activity  such  as  fly  ash.  The  behavior  of  soluble  aerosols  will  be  treated  ' 
in  Section  6  on  relative  humidity  effects. 


3  DESCRIPTIVE  CASE  STUDIES 

In  this  section  we  present  descriptive  aerosol  data  from  three 
research  programs.  These  data  were  selected  because  they  describe  many  of 
the  aerosol  properties  pertinent  to  infrared  extinction,  provide  information 
on  temporal  and  spatial  variability,  and  describe  certain  processes  that 
influence  the  space/time  distribution  of  aerosol  properties.  Case  I  pre¬ 
sents  a  global  view  of  aerosol  distributions  and  illustrates  how  aerosol 
properties  differ  among  various  airmass  types.  Case  II  presents  data  on  the 
marine  aerosol  and  includes  a  study  of  the  behavior  of  windblown  dust. 

Case  III  describes  the  behavior  and  characteristics  of  certain  urba»  aerosols. 

Case  I  Global  Monitoring  Station  Data 

Bigg  (1980)  published  a  comparison  of  the  aerosol  characteristics 
at  four  global  baseline  atmospheric  monitoring  stations.  The  stations  are 
to  establish  reference  levels  of  atmospheric  constituents  that  may  affect 
global  climate.  The  data  from  these  stations  were  found  useful  for  our 
purpose  because  features  of  the  aerosol  at  each  station  and  the  differences 
among  them  illustrate  many  pertinent  features  of  aerosol  variability. 

The  four  stations  were  Point  Barrow,  Alaska;  Mauna  Loa  Observatory, 
Hawaii;  Cape  Grim,  Tasmania;  and  a  South  Pole  station.  Table  3  provides 
descriptive  site  data.  Periods  of  observation  at  each  site  are  indicated, 
along  with  a  brief  description  of  certain  pertinent  meteorological  conditions. 

Fig.  5  shows  representative  number  density  versus  particle  size  for 
the  aerosol  at  each  of  the  four  sites.  Fig*  6  is  a  similar  plot  of  the  mass 
distribution.  Note  the  following  features  of  Fig.  S. 

(1)  Number  densities  in  the  size  range  from  0.03  pm  to  0.08  pm  are 
quite  similar  among  the  four  stations.  Bigg  reported  that  the  size  distribution 
of  the  aerosol  varied  little  with  time  at  Point  Barrow.  He  found  that  the 
predominant  chemical  constituents  in  the  small  particle  component  at  Point  Barrpw 
were  sulferic  acid  and  ammonium  sulfate.  He  attributed  the  relatively  high 
number  concentrations  above  0.2  ym  to  a  sea-salt  component.  He  postulated 
that  the  small  particles  were  transported  poleward  from  mid-latitude  regions 
by  the  high-altitude  circulation.  The  uniform  number  densities  of  the  small 
particles  and  the  remote  origin  of  the  small  particles  at  Point  Barrow  demon¬ 
strate  the  permanence  and  long  residence  time  of  the  small  particle 
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Fig.  5.  The  number  density  distributions  of  particles  at  (a)  Barrow, 
Alaska,  (b)  Mauna  Loa,  Hawaii,  (c)  Cape  Grim,  Tasmania  and 
(d)  South  Pole  (Bigg,  1980} . 


component  and  indicate  a  more  or  less  steady-state  global  production  mechan¬ 
ism.  Sulfur  compounds  dominated  the  South  Pole  aerosol  with  a  much  smaller 
sea-salt  component.  Bigg  attributes  the  origin  of  the  South  Pole  aerosol  to 
mechanisms  similar  to  those  ascribed  to  the  Point  Barrow  aerosol.  The  low 
concentration  of  particles  aboye  0.2  pm  at  the  South  Pole  may  be  due  to 
washout  by  precipitation  or  to  the  high  elevation  of  the  site. 

(2)  Particle  concentrations  in  the  Mauna  Loa  aerosol  were  low  at  all 
particle  sizes  above  0.1  pm.  '’his  result  differs  perhaps  from  the  type  of 
aerosol  expected  in  a  maritime  location.  Absence  of  the  large  particles  is 
due  to  the  high  elevation  of  the  Mauna  Loa  site  and  the  inhibiting  effects  of 
the  trade  wind  inversion  on  vertical  transport  of  large  particles.  The  Mauna 
Loa  aerosol  features  are  characteristic  of  aerosols  often  found  in  the  tropo¬ 
sphere  above  the  boundary  layer.  Junge  (1972)  reported  similar  features  in 
aerosols  above  the  boundary  layer  in  the  North  Atlantic.  The  chemical  compo¬ 
sition  at  Mauna  Loa  was  similar  to  that  of  the  high  latitude  aerosols,  with 
a  predominance  of  water  soluble  sulfur  compounds. 
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(3)  The  distinctive  feature  of  the  Cape  Grim  aerosol  is  the  high 
particle  count  in  the  size  range  above  0.2  ym.  Mote  tliat  in  Bigg's  study, 
the  upper  size  limit  of  distributions  reported  is  about  1  pm.  Since  the 
Cape  Grim  aerosol  is  clearly  a  marine  aerosol,  it  is  reasonable  to  assume 
that  signifi.  ant  number  concentrations  would  be  found  for  particles  as  large 
as  20  ym.  The  small  particle  component  of  the  Cape  Grim  aerosol  was  dominated 
hy  sulfur  compounds.  The  large  particle  component  was  almost  entirely  sea 
salt. 

The  particle  mass  distribution  plotted  in  Fig.  6  shows 
several  important  features.  Note  the  concentration  of  mass  in  very  small  par¬ 
ticles  in  the  Mauna  Loa  aerosol.  The  mass  is  distributed  among  larger  parti¬ 
cles  ranging  from  .04  ym  to  0.5  ym  for  the  Point  Barrow  and  South  Pole  aero¬ 
sols.  The  Cape  Grim  aerosol  shows  a  large  mass  distribution  at  much  larger 


Fig.  6.  The  mass  distributions  of  particles  at  (a)  Barrow,  Alaska,  (b) 
Mauna  Loa,  Hawaii,  (c)  Cape  Grim,  Tasmania  and  (d.  South  Pole 
(Bigg,  1980) . 
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particle  sizes.  A  significant  mass  concentration  at  Cape  Grim  is  likely 
above  the  measurement  size  limit  at  1  pm.  Bigg  estimated  the  total  mass 
concentration  (integrated  over  particle  sizes  to  1  pm)  as  follows: 

Mauna  Loa 
Point  Barrow 
South  Pole 
Cape  Grim 

The  particle  size  distributions  and  mass  concentrations  at  three  of 
these  stations  clearly  permit  classification. of  their  aerosols  according  to 
the  LOWTRAN  models  as  follows: 

Mauna  Loa  Tropospheric 

Point  Barrow  Rural 

Cape  Grim  Maritime 

The  South  Pole  aerosol  falls  between  the  rural  and  tropospheric  models 

Another  important  feature  of  Fig.  6  is  the  modal  distribution  of  the 
aerosol  mass.  The  Mauna  Loa  aerosol  shows  two  modes,  but  both  are  in  the 
small -particle  component .  The  Point  Barrow  and  South  Pole  aerosols  show  a 
weak  secondary  mode  near  0.4  pm.  The  Cape  Grim  aerosol  clearly  shows  both 
the  accumulation  mode  and  the  coarse-particle  mode.  In  this  case,  the  small- 
particle  mode  represents  the  stable  aerosol  with  a  continental  origin,  while 
the  coarse-particle  mode  represents  the  sea-salt  aerosol. 

Figure  7  shows  electron  microscope  photographs  of  aerosol  samples 
taken  at  the  four  sites.  These  photographs  clearly  show  the  different  size 
distributions  among  the  aerosols. 

Case  II  North  Atlantic  Data  from  the  R.  V.  Meteor 

Junge  and  Jaenicke  (1971)  and  Junge  (1972)  contain  reports  of  the 
voyage  of  the  German  ship  R.  V.  Meteor  in  the  North  Atlantic  in  1969.  Fig. 

8  is  a  map  of  the  expedition.  Extensive  aerosol  measurements  were  taken  on 
the  northbound  leg  of  the  voyage  from  April  14  to  May  7.  The  ship  crossed  the 
Intertropical  Convergence  zone  (XTCZ)  on  the  northbound  leg  on  April  17th. 

The  airmass  south  of  the  ITCZ  was  considered  maritime  in  character.  The 
airmass  immediately  north  of  the  ITCZ  had  its  origin  in  the  Sahara  desert 
and  was  laden  with  Sahara  dust. 


.08  pg  ,m 
0.8  pg  m  ■* 
0.14  pg  m  ^ 
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Pig.  7.  Particles  collected  with  1  sun  jet  impactor  at  (a)  Barrow  Alaska,  1 

(b)  Mauna  Loa,  Hawaii,  (c)  Cape  Grim,  Tasmania  and  (d)  South  | 

Pole.  The  fields  of  view  shown  correspond  roughly  to  100  cnr 
of  sampled  air.  The  inset  is  2  urn  in  length  (Bigg,  1980) . 

■  ’  • 

The  visual  range  was  recorded  at  the  green  wavelength  (0.55  pm)  with  a  | 

specially  designed  instrument.  Traces  from  this  device  indicated  a  reduction  in 
visual  range  from  100  km  to  20  km  upon  crossing  the  ITCZ.  The  authors  reported  \ 

that  dust  was  deposited  on  the  ship’s  deck  after  crossing  the  ITCZ.  Further  i 

north,  the  airmass  returned  to  a  purely  maritime  character  on  about  April  24. 

Fig.  9  shows  a  twice  daily  analysis  of  the  trajectory  of  the  a* -mass  observed 

ait  each  point  on  the  northbound  leg  after  crossing  the  ITCZ.  Successive  ! 

dpts  on  each  trajectory  represent  estimated  12  hour  airmass  movement. 

Fig.  10  shows  time  traces  of  the  output  from  various  instruments 
aboard  ship.  The  instruments  are  named  and  their  output  is  described  in  the 
caption.  The  following  features  of  the  traces  provide  ins  light  into  the 
variation  of  several  aerosol  properties. 
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Fig.  8.  Map  of  the  Atlantic  expedition  of  the  German  Research  Vessel 

Meteor,  Spring  1969.  Dates  are  indicated  at  various  positions 
in  the  expedition  by  day. month  (Junge  and  Jaenicke,  1971) . 

Trace  A  shows  total  couhts  of  the  Aitken  nuclei  in  the  raijge  .001  ym 
to  .01  ym.  Note  that  the  number  concentration  varies  little  with  airmass 
changes,  with  a  mean  value  about  600  cm  The  authors  reported  surprise  at 
the  steady  values  and  their  unexpectedly  high  magnitudes.  These  data  were 
the  first  real  evidence  of  a  steady  state  production  mechanism  to  offset  the 
loss  of  particles  in  this  size  range  through  coagulation.  The  authors 
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Legend 


A  :  Aitk.cn  condensation  nuclei 
R  :  Royco  particle  counter  (radius  >  0.3  pm) 
I  i  Single  stage  impactor  (radius  >  3  pm) 

PI  :  Rotating  Impactor  (radius  >  19  pm) 

S  :  Visibility  with  green  wavelength 
<  :  Wind  speed 


For  I  and  FI 

_  Soluble  particles 

^  Insoluble  particles 
Sum  of.  both 


Pig.  10.  Particle  concentrations  measured  with  various  instruments  as  a 

function  of  time  (GMT)  and  position  of  the  ship  on  30°  W.  Hatched 
areas  indicate  wind  direction  80°  -  120°  relative  to  main  ship  axis 
with  high  probability  of  pollution  from  the  ship  (Junge  and  Jaenickc 
197 


considered  this  nuclei  count  as  low,  however,  compared  to  those  found  in 
continental  regions. 

Trace  R  represents  total  number  concentrations  in  the  range  of  par¬ 
ticle  sizes  above  0.3  yzn.  Note  the  sharp  increase  in  number  concentrations 
on  April  17  upon  moving  into  the  Saharan  airmass.  This  increase  in  number 
concentration  is  clearly  the  cause  of  the  simultaneous  reduction  in  visi¬ 
bility  as  shown  in  Trace  S.  Number  concentrations  in  this  size  range  re¬ 
turned,  to  much  lower  levels  on  April  24  when  the  ship  again  entered  an  air- 
mass  with  a  long  maritime  history.  The  subsequent  mean  value  of  visibility 
was  significantly  higher  than  that  observed  in  the  Saharan  airmass  but  was 
also  quite  variable.  The  variation  shows  a  strong  inverse  correlation  with 
the  number  concentrations  in  Trace  R. 

Trace  Z  represents  total  number  concentrations  for  particle  sizes 
greater  than  3  pm.  This  trace  shows  the  counts  separately  for  the  soluble 
and  the  insoluble  components.  Trace  I  shows  a  sharp  increase  in  the  counts 
upon  entry  into  the  Saharan  airmass  on  April  17,  witn  most  of  the  increase 
in  the  insoluble  component  (dashed  line) .  These  are  the  Sahara  dust  parti¬ 
cles.  Upon  entry  into  the  maritime  air  on  April  24  the  total  particle  count 
in  the  range  above  3  pm  decreased  markedly.  The  counts  for  the  insoluble 
component  decreased  by  almost  a  factor  of  100.  The  soluble  component  domi¬ 
nated  the  distribution  after  April  25  with  counts  about  the  same  as  those 
found  for  the  soluble  component  before  that  date. 

Trace  FI  shows  that  the  number  concentrations  of  particles  above  19  pm 
are  three  to  four  orders  of  magnitude  lower  than  the  counts  of  particles  above 
3  pm.  Because  of  their  low  number  density,  their  existence  in  the  undis¬ 
turbed  marine  aerosol  is  of  no  consequence  to  the  strength  of  infrared 
extinction.  Me  refer  the  reader  to  the  source  paper  for  a  discussion  of  the 
properties  of  these  very  large  aerosols. 

Junge  reported  a  detailed  analysis,  based  in  part  on  the  R.  V.  Meteor 
data,  of  the  behavior  of  sea-salt  aerosols  in  the  marine  environment.  The 
following  are  pertinent  conclusions  of  his  analysis: 

(1)  The  upper  limit  of  the  size  distribution  is  20  pm  in  the 
sea-salt  aerosol.  He  concluded  that  particles  larger  than  20  pm  would  not 
be  found  in  significant  numbers  sxcept  in  the  layer  (on  the  order  of  10  meters 
thickness)  immediately  above  the  ocean  surface.  The  particles  larger  than 
20  pm  are  produced  by  sea  spray  and  are  lost  rapidly  by. sedimentation. 
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(2)  The  lower  size  limit  of  sea-salt  particles  is  about  0.1  yin. 
Although  their  number  concentrations  may  be  quite  high  in  the  size  range 
from  0.2  ym  to  1.0  ym,  their  counts  cure  an  order  of  magnitude  less  than  those 
of  other  soluble  compounds  in  the  background  aerosol  in  this  size  range. 

(3)  The  chemical  composition  of  sea -salt  particles  is  dominated  by 
sodium  chloride.  However,  significant  amounts  of  other  soluble  and  insoluble 
compounds  are  found  within  the  sea-salt  crystals. 

(4)  Particle  counts  decrease  rapidly  with  height  above  about  0.5  km 
over  the  ocean,  with  practically  no  sea  salt'  found  above  2-3  km.  This  finding 
does  not  necessarily  hold  over  land  because  vertical  transport  mechanisms 
operate  more  effectively  over  land. 

(5)  Residence  times  for  sea-salt  particles  in  the  size  range  below 
10  ym  average  about  three  days.  The  particles  above  10  ym  will  have  shorter 
residence  times.  Residence  time  depends  both  on  sedimentation  and  on  rain- 
out  and  washout. 

Junge’s  work  represents  one  of  the  early  efforts  to  understand  the 
global  background  aerosol.  He  found  that  the  composition  of  the  background 
aerosol,  except  for  the  sea-salt  component,'  was  very  uniform  in  both  the 
horizontal  and  the  vertical  in  oceanic  regions.  He  reached  the  tentative 
conclusion  that  sulfate  compounds  made  up  more  them  50  percent  of  these 
aerosols.  The  data  presented  by  Bigg  in  Case  I  are  results  of  just  one  of 
several  later  research  programs  that  have  confirmed  his  findings. 

Junge  found  that  the  principal  constituent  of  the  global  background 
aerosol  is  ammonium  sulphate  These  particles  are  known  to  have  a  continental 
origin.  Their  number  concentrations  dominate  the  particle  distribution  in  the 
size  range  from  0.05  ym  to  0.5  ym.  Fig. -11  shows  the  relationships  between 
the  continental  and  marine  aerosols  in  this  size  range. 

Case  III  Urban  Aerosols 

Willeke,  et  al.  (1974)  and  Willeke  and  Whitby  (1975)  are  reports  of 
extensive  studies  of  urban  aerosols  in  the  Denver,  CO  area,  as  well  as  mea¬ 
surements  at  a  number  of  other  locations.  Measurements  were  taken  at  two 
different  locations  in  the  Denver  area  in  the  fall  of  1971.  Willeke  and 
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Fig.  11.  Idealized  size  distributions  of  continental  and  marine  aerosols. 

(Adapted  from  Junge,  1972). 

Whitby  compared  these  measurements  to  those  obtained  at  other  locations  dur¬ 
ing  the  period  from  1969-74.  All  measurements  were  taken  with  the  Minnesota 
Aerosol  Analyzing  System  (MAAS) .'  Their  results  represent  several  thousand 
in-situ  measurements  of  aerosol  particle  size  distributions.  Findings  perti¬ 
nent  to  the  problem  of  visible  and  infrared  extinction  are  presented  here. 

The  MAAS  measures  particles  in  the  size  range  from  .007$  pm  to 
17.5  pm  diameter.  Surface  area  distributions  were  computed  by  assuming  a 
spherical  shape  for  the  particles.  Visual  range  was  determined  from  nephelo- 
meter  measurements  of  the  scattering  coefficient  in  the  visible  spectrum. 

Fig.  12  shows  aerosol  particle  surface  area  distributions  as  a  func¬ 
tion  of  particle  size  under  five  different  conditions,  fable  4  lists  a  number 
of  key  parameters  for  these  five  cases. 

Variation  in  the  visual  range  in  the  various  types  of  aero¬ 
sols  is  a  matter  of  interest  here.  In  the  five  aerosol  data  sets  presented 
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Moderate  urban  pollution  plus  heavy  traffic.  Denver ,  CO;  2*s 
km  West  of  downtown,  Oct.  26,  1971,  1640-1700  hours. 

Heavy  urban  pollution.  Denver-Welby,  CO;  10km  NE  of  downtown; 
Nov.  9,  1971,  0700-0900  hours. 

Heavy  urban  pollution.  Denver-Welby,  CO;  10  km  NE  of  downtown; 
Nov.  9,  1971 ,  0920-1330  hours. 

Urban  district  background  plus  anthropogenic  Denver-Welby,  CO; 
10  km  NE  of  downtown;  0200-  0400  hours. 

Average  western  continental  background.  Port  Collins,  CO; 

Aug.  13-14,  1970. 


Pig.  12.  Aerosol  surface  area  distributions  under  the  five  conditions 
described  in  the  legend.  This  figure  was  adapted  from  data  presented  by 
Willeke  and  Whitby.  The  plots  for  distributions  4  and  5  are  indistin¬ 
guishable  beyond  a  diameter  of  0.5  jim  in  this  presentation  (Willeke  arid 
Whitby,  1975). 
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in  Fig.  12,  only  the  first  included  a  report  of  visual  range.  In  order  to 
obtain  comparative  data  on  visual  range,  additional  cases  reported  by  Willeke 
and  Whitby  were  analyzed  in  search  of  comparable  size  distributions  accom¬ 
panied  by  visual  range  data.  The  right  hand  side  of  Table  4  contains  data 
for  three  cases  containing  reports  of  visual  range.  The  size  distribution 
in  each  of  these  three  cases  was  very  similar  to  that  of  the  corresponding 
case(s)  on  the  left  side  of  the  table. 

Note  from  Fig.  12  that  in  most  urban  aerosols,  the  particle  surface 
area  distributions  are  several  orders  of  magnitude  higher,  in  some  subrange 
below  0.5  pm,  than  corresponding  values  in  the  continental  background  aero¬ 
sols.  The  specific  subrange  in  which  this  occurs  depends  on  the  nature, 
source  and  age  of  the  pollutants  and  has  a  direct  bearing  on  the  effect  of 
the  urban  aerosol  on  visual  range. 

We  note  that  in  Case  1,  the  visual  range  is  not  appreciably  less 
than  in  Cara  4,  in  spite  of  the  very  high  particle  surface  area  in  the 
subrange  below  0.1  pm.  The  small  particles  in  this  subrange  are  not  effi¬ 
cient  attenuators  of  visible  light.  Cases  2  and  3  show  a  visual  range  much 
lower  (deduced  from  the  Pomona,  CA  data)  than  that  reported  in  Case  1.  We 
will  show  in  Section  5  why  the  small  particles  in  Case  1  have  little  effect 
on  the  visual  range. 

Note  that  the  size  distributions  in  Fig.  12  show  that  Cases  1  and  4 
have  similar  number  concentrations  in  the  size  range  from  0.1  pm  to  1.0  pm. 
The  visual  range  values  are  similar  in  these  two  cases.  In  fact,  the  similar 
values  of  visual  range  reflect  the  similarity  of  particle  counts  in  the  size 
range  that  controls  visibility  restriction.  The  numerous  particles  in  the 
nuclei  mode  in  Case  1  had  a  small  effect  on  visibility. 

The  downtown  Denver  site  is  located  about  100  meters  from  a  freeway. 
Note  that  the  period  for  Case  1  was  during  evening  rush  hour  traffic.  The 
observation  site  was  downwind  from  the  freeway  in  a  brisk  30  km  hr”^  wind. 

The  particle  size  distribution  below  0.1  pm  represents  fresh  combustion 
aerosols  generated  in  the  freeway  traffic.  Vehicular  traffic  in  the  City 
Maintenance  Yard  had  subsided  for  the  day,  and  few  power  plant  emissions 
reached  the  site;  hence  the  particle  counts  in  the  range  above  1  pm  were 
very  low. 
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The  high  value  of  surface  area  distribution  below  0.1  pm  in  Case  1 
is  a  transient  feature.  Within  one  or  two  hours,  the  Aitken  nuclei  coagulated 
and  formed  another  peak  in  the  distribution  around  0.1  pm.  Case  2  illustrates 
the  distribution  measured  on  another  date  at  a  site  located  some  6-7  km  from 
the  nearest  freeway.  This  distribution  is  similar  to  that  expected  from 
Case  1  after  the  coagulation  process  has  acted  for  a  few  i;ours.  Case  3  repre¬ 
sents  further  aging  of  a  similar  aerosol,  just  two  hours  after  Case  2  at  the 
same  location.  These  temporal  patterns  are  characteristic  of  heavy  urban 
pollution. 

The  coarse  particles  at  10  pm  in  Cases  2  and  3  are  unrelated  to  the 
fine  particles.  These  particles  were  attributed  to  settling  out  of  fly  ash 
from  power  plants. 

Case  4  is  a  surface  distribution  characteristic  of  the  urban  district 
background.  This  represents  a  condition  of  aged  aerosols  in  which  dispersion 
has  decreased  the  number  concentrations.  This  distribution  shows  two  peaks. 
The  peak  at  0.08  pm  represents  the  coagulation  of  fresh  nuclei,  while  the 
second  peak  represents  the  remains  of  older  aerosols. 
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4  PROCESSES  AFFECTING  THE  NATURE  OF  THE  AEROSOL 


The  total  mass  content  of  aerosol  in  the  atmosphere  and  its  varia¬ 
tion  in  time  is  determined  by  production  and  loss.  The  spatial  distribu¬ 
tion  of  the  aerosol,  its  chemical  composition,  and  its  size  distributions 
are  determined  by  various  processes.  This  section  discusses  some  processes 
of  importance  to  the  problem  of  selecting  an  appropriate  aerosol  model. 

Fig.  13  summarizes  the  definitions  of  the  particle  size  ranges  for 
the  three  modes  commonly  found  in  aerosol  size  distributions.  It  identifies 
the  principal  aerosol  sources  and  characterizes  the  residence  times  typical 
of  the  aerosols  in  each  mode. 

Fin*  particles  -4*.  Coaru  particle* 
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Fig.  13.  Nomenclature,  origin,  lifetimes,  and  sizes  of  aerosols  that  con¬ 
tribute  to  the  three  modes  observed  in  the  atmosphere.  Size  ranges  refer 
to  particle  diameter.  (Willeke  and  Whitby,  1975) 

Production  Processes 

We  may  classify  aerosol  production  mechanisms  according  to  the  time 
scale  on  which  they  influence  the  atmospheric  aerosol.  The  time  scale  of 
interest  to  us  lies  in  the  range  of  a  few  hours  to  a  few  days.  Many  im¬ 
portant  production  processes  affect  the  aerosol  on  a  much  longer  time 
scale.  Their  products  form  the  global  background  aerosol  with  rather 
uniform  spatial  and  temporal  characteristics.  Prom  our  perspective  these 
processes  act  as  more  or  less  steady  state  sources  for  maintenance  of  the 
background  aerosol.  Their  products  sure  represented  in  LOWTRAN  by  the  rural 
aerosol  model  in  the  boundary  layer  and  by  the  tropospheric  aerosol  model 
in  the  atmosphere  above  the  boundary  layer.  The  following  are  the  prin¬ 
cipal  processes  that  maintain  the  background  aerosol. 


-  Fossil  fuel  combustion  and  forest  fires  produce  global  background 
aerosols  principally  in  the  nuclei  range.  Certain  aerosols  in  the  accumula¬ 
tion  range  are  also  direct  products  of  combustion.  The  aerosols  in  the 
nuclei  range  coagulate  rapidly  to  form  aerosols  in  the  accumulation  range. 
Both  direct  and  indirect  products  have  long  atmospheric  residence  times. 
These  processes  are  the  origin  of  many  of  the  sulfur,  nitrogen  and  carbon  , 
compounds  in  the  background  aerosol. 

-  Plant  emissions  produce  aerosols  in  both  the  nuclei  and  accumula¬ 
tion  ranges.  Certain  organic  aerosols  and  nitrogen  compounds  result  from 
this  source.  Hie  product  of  these  emissions  is  often  seen  as  blue  haze  in 
regions  of  extensive  vegetation. 

-  Sea  spray  and  sea  foam  produce  aerosols  across  a  wide  size  range. 
Those  in  the  range  below  0. 3-0.5  pm  become  a  part  of  the  background  aerosol. 
Their  number  concentrations  in  this  size  range  are  much  lower,  however, 
than  those  of  aerosols  from  other  sources. 

Certain  production  processes  result  in  significant  variability 
in  aerosol  properties  on  the  tiro  scale  of  hours  to  days.  The  following 
are  the  principal  mechanisms  of  interest  to  our  forecasting  problem. 

-  Fossil  fuel  combustion  in  urban  and  industrial  complexes  pro¬ 
duces  short  term  local  variability.  Aerosols  result  in  all  three  size 
ranges.  Those  in  the  nuclei  and  accumulation  ranges  ultimately  enter  the 
background  aerosol.  Large  particle  emissions,  such  as  fly  ash,  are 
strongly  influenced  by  sedimentation  and  by  precipitation  scavenging.  As 
a  result,  they  exhibit  short  residence  times. 

-  Sea  spray  and  sea  foam  are  the  primary  production  mechanism  for 
sea-salt  particles.  Those  in  the  size  range  from  2-5  pm  exhibit  residence 
times  of  a  few  days.  The  larger  sea-salt  particles  (above  10  pm  are  removed 
by  sedimentation  on  the  time  scale  of  a  few  hours.  Fig.  14  is  taken  from 
data  published  by  Wells  et  al  (1977).  They  show  the  dependence  of  coarse 
particle  concentrations  in  maritime  air  on  wind  speed.  The  very  large 
particles  (above  20  pm)  in  these  distributions  are  found  only  in  a  shallow 
layer  (on  the  order  of  10  meters  in  depth)  above  the  ocean  surface.  They 
remain  airborne  for  time  periods  of  less  than  a  few  hours  because  of 
sedimentation.  The  depth  of  the  marine  aerosol  layer  (for  particles 
smaller  than  20  pm)  depends  on  the  strength  of  vertical  turbulent  mixing. 
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Fig.  14.  Characteristic  aerosol  number  density  distributions  in  maritime 
air  as  a  function  of  velocity  for  a  RH  of  80  percent  (solid 
lines) .  Predictions  of  the  maritime  component  of  the  Hells 
model  based  on  these  data  are  also  shown  (dotted  lines) (Hells, 
1977) . 

The  sea-salt  aerosol  is  of  primary  importance  in  forecasting' infrared 
extinction.  The  size  of  the  sea-salt  particles  and  their  response  to 
changes  in  relative  humidity  are  the  largest  overall  source  of  time/ 
space  variability  in  the  large  particle  concentrations.  An  estimate 
of  their  relative  concentration  is  quite  important  in  prediction  of 
aerosol  extinction.  The  behavior  and  importance  of  the  sea-salt  aerosol 
will  be  discussed  more  in  Sections  5  and  6. 


Transport  and  Diffusion 


Transport  and  diffusion  are  both  effective  mechanisms  in  redistri¬ 
buting  the  aerosol  mass.  Long-range  transport  of  the  small  particles  in 
the  accumulation  mode  accounts  for  the  globally  uniform  background  aerosol. 
Transport  also  plays  an  important  role  on  the  time  scale  of  hours  ;o  days 
in  carrying  the  sea-salt  in  maritime  airmasses  into  continental  regions. 
The  effect  of  transport  of  the  aerosol  into  a  region  is  counterbalanced  in 
both  cases  by  mechanisms  that  remove  the  aerosol  from  the  atmosphere. 
Transport  is  important  over  short  time/distance  scales  in  relocation  of 
urban  aerosols  downwind  from  their  source. 

Diffusion  plays  an  important  small-scale  role  in  the  redistribu- 
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tion  of  urban  aerosols.  The  general  effect  of  diffusion  is  to  decrease  their 
number  concentrations.  Vertical  diffusion  acts  to  carry  urban  aerosols 
into  the  free  atmosphere  for  rapid  dispersion.  Horizontal  diffusion  acts 
to  reduce  their  number  concentrations  while  it  increases  their  area  of 
influence. 

Removal  Mechanisms 

Threw  principal  mechanisms  for  removal  of  aerosols  are  sedimentation, 
rainout  and  washout.  Sedimentation  acts  to  a  small  degree  on  the  accumula¬ 
tion  mr.de .  However,  its  greatest  importance  is  in  the  removal  of  large 
particles. 

Kasten  (1968)  presented  equations  for  falling  speed  of  aerosol 
particles.  His  development  considered  only  gravitational,  viscoiis,  and 
buoyant  forces.  From  his  equations,  the  terminal  falling  speed  of  an 
aerosol  particle  in  the  size  range  from  1  ym  to  10  ym  is  given  by  the 
relation: 

w«appr2. 

where  w  is  the  falling  speed,  a  is  a  proportionality  constant,  p  is  the 

P 

density  of  the  particle  and  r  is  particle  radius.  This  relation  shows 
that  sedimentation  rapidly  becomes  more  effective  with  increasing  particle 
size,  and  falling  speed  is  directly  proportional  to  particle  density. 

Junge  (1972)  performed  an  extensive  analysis  of  expected  residence 
times  of  sea-salt  aerosols  over  the  ocean.  He  based  his  study,  in  part, 
on  data  by  Eriksson  (1957) ,  who  studied  the  balance  between  sea-salt 
aerosol  production  and  loss  by  sedimentation.  Data  are  presented  in  Table 
5,  with  expected  residence  time  as  a  function  of  particle  mass. '  We  have 
converted  particle  mass  to  equivalent  diameter  for  three  different  particle 
densities. 

Eriksson's  first  estimate  is  based  on  gravitational  fall 
velocities  alone.  The  data  show  very  long  residence  times  for  particles 
in  the  range  1-2  ym.  Eriksson  revised  this  estimate  to  reconcile  the 
estimated  residence  time  with  sea-salt  production  estimates  and 
observed  sea-salt  concentrations.  Junge  concluded  from 


Erihss on's  estimates  that  sedimentation  could  not  be  the  primary  removal 
mechanism  for  particles  smaller  them  10  pm  diameter.  He  did  conclude  that 
for  particles  larger  than  10  pm,  sedimentation  acting  alone  would  result 
in  residence  times  considerably  less  than  one  day.  He  attributed  the  short 
residence  times  observed  for  the  smaller  particles  to  removal  by  rainout 
and  washout. 


Table  5.  Estimates  of  residence  times  for  sea-salt  particles  as  a  function 
of  particle  mass.  The  equivalent  diameters  are  my  values  cal¬ 
culate^  by  assuming  spherical  particles  wijjh  a  density  of  2.1 
gm  cm  (dry  sodium  chloride),  1.3  gm  cm  (sodium  chloride  in 
solution,  and  1.0  gm  cm  (pure  water).  (Junge,  1372,  Eriksson, 
1957) 


Mass  (grams) 

io-12 

io-11 

io'10 

io*9 

Residence 

time  (days) 

Eriksson's  estimate 

82 

16 

0.6 

0.5 

(sedimentation  alone) 

Eriksson's  estimate 

3.5 

1.0 

0.6 

0.5 

(production  and 
sedimentation ) 

Junge ’s  estimate 

1.9 

1.6 

1.0 

.26 

(see  text) 

' 

Equivalent  diameter  (pm) 

■HBHB 

■ 

Density  *  2.1  gm  cm  3 

BSNFfl 

3.0 

6.5 

14 

Density  =*  1.3  gm  cm-3 

3.6 

7.7 

17 

Density  =  1.0 

3.9 

8.4  . 

18 

Rainout  occurs  when  aerosol  particles  act  as  cloud  condensation 
nuclei  and  subsequently  fall  out  in  precipitation.  Washout  occurs  when 
falling  precipitation  particles  pick  up  aerosol  particles  during  their 


fall.  Precipitation  scavenging  usually  refers  to  washout;  however,  in 
field  measurements  of  precipitation  scavenging,  it  is  often  difficult  to 
separate  the  effects  of  rainout  from  those  of  washout. 

Several  theoretical  studies  have  shown  the  theoretical  washout 
efficiency  of  precipitation.  This  efficiency  is  a  complex  function  of  both 
the  aerosol  size  distribution  and  the  raina,  op  size  distribution.  We 
present  one  set  of  theoretical  results  to  illustrate  the  important  features 

Dana  and  Hales  (1976)  studied  the  theoretical  dependence  of  wash¬ 
out  efficiencies  on  both  aerosol  and  raindrop  size  distributions.  Fig.  15 
shows  the  principal  features  pertinent  to  radiation  extinction.  Note 
that  both  axes  are  a  logarithmic  scale.  Collection  efficiency  curves 
are  shown  for  two  raindrop  size  distribution  in  which  it  was  assumed  that 
all  raindrops  were  the  same  size.  Fig.  15  contains  three  regions  in  the 
aerosol  size  spectrum  that  represent  three  different  aerosol  collection 
processes,  tn  the  region  below  0.05  urn  collection  occurs  by  Brownian 
diffusion.  From  0.05  pm  to  about  1  pm  collection  occurs  by  impaction., 

This  process  consists  of  inertialess  flow  in  which  the  aerosol  particles 
are  assumed  to  be  deflected  by  the  airstream  around  the  raindrop.  In 
the  region  above  1.0  pm,  collection  occurs  in  inertial  flow  in  which  it 
is  assumed  that  aerosol  particle  position  is  unaffected  by  the  airstream 
around  the  raindrop.  We  find  three  noteworthy  features  in  Fig.  15:  (1) 
very  high  collection  efficiencies  above  1  pm;  (2)  minimum  in  collection 
efficiency  around  0.05  pm.  (3)  Collection  efficiences  much  higher  (about 
two  ordersof  magnitude)  for  the  large  particles  than  for  the  small  parti¬ 
cles. 

RadJce  et  al.  (1980)  report  the  results  of  measurements  of  aerosol 

washout  in  various  locations. _ Table  6  lists  descriptive  data  on  the 

measurement  program.  Fig.  16  and  17  show  aerosol  size  distributions  at 
two  sites  before  and  after  a  rain  shower  event.  Note  the  well  defined 
minimum  in  the  washout  around  0.5-0. 8  pm  particles  diameter  in  both  cases. 

Fig.  18  shows  collection  efficiency  as  a  function  of  aerosol 
particle  size  for  each  event  listed  in  Table  6.  Note  that  the  vertical 
coordinate  in  Fig .  18  represents  the  percentage  of  particles  removed  by 
scavenging  while  that  in  Fig.  15  represents  a  parameter  named  the  SC2'  ng- 
ing  collection  efficiency. 


RADIUS  (ym> 

Fig.  15.  Theoretical  collection  efficiencies  (assuming  complete  retention) 
used  in  calculations,  and  selected  experimental  results.  R 
refers  to  the  raindrop  size.  Experimental t  x.  Starr  and  Mason 
(R  -  0.05  cm)  (1966) :  0.  Walton  and  Woodcock  (R  -  0.05  cm)  (1970) : 

+.  Adam  and  Semonin  (R  in  ())(1970):  A.  Sood  and  Jackson  (R  -  0.05 
cm) (1970) :  A.  Engelmann  (R  -  0.02  cm) (1965).  The  abscissa  is 
aerosol  particle  radius.  (Adapted  from  Dana  and  Hales,  1976) . 
Sources  of  experimental  data  are  identified  in  the  source  paper. 


Table  6.  Conditions  under  which  aerosol  precipitation  scavenging  measurements 
were  obtained  (Radke  et  al.,  1980). 
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These  parameters  are  not  the  same;  we  refer  the  reader  to  the  source  papers 
for  a  definition  of  the  distinction.  Both  are  indicators  of  washout  effi¬ 
ciency,  however. 

The  following  features  in  Fig.  16,  17  and  18  are  noted. 

(1)  A  pronounced  minimum  in  washout  efficiency  appears  in  every 
case.  This  minimum  appears  in  the  range  of  aerosol  sizes  from  0.4  to  1.1 
pm  diameter.  Compare  this  to  the  theoretical  minimum  in  Fig.  15  at  0.05 
pm  radius . 

(2)  Scavenging  efficiency  is  high  in  every  case  for  large  parti¬ 
cles.  Effectively,  it  reaches  80  percent  at  2-4  pm  particle  diameter  and 
100  percent  for  particles  larger  than  10-12  pm.  This  result  agrees  well 
with  theoretical  results  shown  in  Fig.  15. 

(3)  The  percentage  of  aerosol  particles  removed  increases  with 
increased  scavenging  time  (May  10,  1976)  for  large  particles  and  decreases 
for  small  particles.  There  is  also  a  pronounced  shift  in  the  scavenging 
gap  toward  small  aerosols  with  increased  scavenging  time. 

(4)  Measured  scavenging  efficiences  for  small  particles  are  much 
greater  than  predicted  from  theory  as  shown  in  Fig.  15.  This  is  believed 
to  result,  at  least  in  part,  from  the  fact  that  many  of  the  small  particles 
act  as  cloud  condensation  nuclei  and  appear  in  the  experimental  data  as 
scavenged  particles.  Theory  of  Fig.  15  does  not  include  this  effect. 

This  phenomenon  may  also  partially  explain  the  shift  in  the  theoretical 
size  for  the  scavenging  gap  toward  larger  particles  in  the  experimental 
data.  Nonetheless,  the  efficiencies  for  small  particles  are  still  sub¬ 
stantially  less  than  those  of  large  particles. 

(5)  The  pronounced  increase  in  scavenging  efficiency  for  large 
particles  begins  at  about  1  pa  to  10  pa  in  both  theoretical  and  experi- 
aental  data. 

Theory  and  experiment  show  that  washout  and/or  rainout  are  very 
important  in  removal  of  aerosol  particles,  particularly  the  large  parti¬ 
cles.  For  a  maritime  aerosol  with  high  particle  number  counts  above  2  pmr 
we  can  expect  that  the  net  effect  6f  precipitation  on  the  size  distribu¬ 
tion  is  to  transform  a  typical  maritime  size  distribution  into  a  dis¬ 
tribution  more  characteristic  of  continental  aerosols. 
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Fig.  16.  Airborne  measurements  of  the  size  spectra  of  aerosol  in  the  plume 

from  the  coal -fired  electric  power  plant  near  Centralia ,  Washington, 
measured  at  5  km  (12  min  travel  time)  downwind  of  the  stack  on  10 
May  1976,  before  and  after  the  plume  was  intercepted  by  a  rain 
shower  (Radke  et  al.,  1980). 
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in  Fig.  16  except  for  "natural"  aerosol  particles  measured, at 
altitude  of  2.5  km  MSL  near. Miles  City,  Montana  on  1  July  1976 
(Radke  et  al. ,  1900) . 
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Pi9*  18.  Percentages  of  aerosol  particles  of  various  sizes  removed  by  pre 
cipitation  scavenging.  See  Table  6  for  details  on  conditions 
under  which  measurements  were  obtained  (Radke  et  al.f  1980). 
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Processes  Affecting  the  Urban  Aerosol 

Dispersion  and  pollutant  source  strengths  are  important  in  the  TDA 
forecasting  problem  in  assessing  the  concentration  of  urban  aerosols. 

Table  7  contains  data  from  selected  sources  on  typical  concentrations  of 
Total  Suspended  Particles  (TSP)  as  well  as  concentrations  of  the  principal 
constituents  in  urban  aerosols.  A  few  reports  also  contain  data  on  the 
fraction  of  the  urban  aerosol  mass  contained  in  the  coarse  particle  range. 
The  size  limit  boundary  between  the  fine  and  coarse  fractions  is  usually 
set  at  either  2.5  urn  or  3.5  pm  in  studies  of  urban  aerosols  because  this 
range  represents  the  upper  size  limit  of  respirable  particles.  The 
following  features  in  Table  7  are  important  for  our  purpose: 

(1)  Total  Suspended  Particulates  (TSP)  (interpreted  as  total  mass 
concentration)  have  values  typically  10-20  times  those  found  in  the  natural 
aerosol. 

(2)  Carbon  compounds  are  an  important  part  of  the  urban  aerosol. 

The  majority  are  organic  compounds.  Some  are  soluble;  some  are  not.  The 
elemental  carbon  (EC)  component  (often  referred  to  as  soot)  is  almost 
entirely  insoluble. 

(3)  Coarse  particles  often  comprise  a  significant  part  of  the 
total  urban  aerosol.  This  component  is  typically  20-60  percent  of  the  total 
mass.  The  coarse  particles  are  made  up  primarily  of  crustal  materials 

(fly  ash  and  soil  derived  aerosols) .  These  materials  are  almost  entirely 
insoluble. 

(4)  He  found  little  data  on  the  rate  at  which  urban  aerosols 
decrease  in  concentration  beyond  the  urban  area  boundary .  Junge  (1963) 
reports  that  concentrations  approach  natural  levels  as  near  as  5  km  from 
the  boundary.  Other  sources  have  indicated  similar  ranges. 

(5)  The  soluble  sulfate  and  nitrate  compounds  are  found  almost 
entirely  in  the  fine  fraction.  Their  behavior  is  similar  to  that  of  the 
background  aerosol;  the  variations  in  total  concentration  are  often 
reflected  in  variations  in  particle  mass  in  that  component. 

In  the  problem  of  LOWTRAN  aerosol  model  selection,  we  based  the 
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Table  7.  Characteristic  properties  of  various  types  of  urban  aerosols  from  selected  sources.  Junge  (1963) 
r^r-orts  a  comparison  the  natural  aerosol.  EC  is  elemental  carbon.  The  boundary  between  fine 

and  coarse  fractions  as  reported  by  Wolff  and  Countess  is  about  3  ym. 


selection  criteria  for  the  urban  model  on  the  assessment  of  the  expected 
concentrations  of  the  urban  aerosol .  This  concentration  at  a  given  time 
will  depend  on  the  strength  of  the  urban  aerosol  sources  and  on  those  mete' 
orological  conditions  which  control  the  effectiveness  of  dispersion  mech¬ 
anisms.  Appendix  A  describes  how  we  arrived  at  selection  criteria  based 
on  a  classification  of  the  total  source  strength  of  an  urban  area  and  on 
boundary  layer  wind  speed  and  stability  conditions. 
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5  FACTORS  THAT  DETERMINE  EXTINCTION  EFFICIENCY 

Aerosol  scattering  at  a  particular  wavelength  is  comnonly  expressed  in 

terms  of  the  scattering  efficiency,  Q  .  The  equation 

5 

1  -  TS  "  «s  nl 

states  the  relationship  for  a  single  particle  of  radius,  r,  in  a  unit  volume. 
Ts  is  transmission;  (1  -  Tg)  represents  the  fractional  depletion,  per  unit  of 
distance  through  the  /olume,  of  the  incident  energy  by  the  particle. 

The  scattering  efficiency  is  related  to  the  scattering  cross- 
section,  o  ,  by 

b 

°s  *  Qs  ** 

The  scattering  cross-section  expresses  the  effectively . scattering 
cross-sectional  area  in  relation  to  the  physical  cross-sectional  area. 

The  scattering  efficiency  defines  the  relationship. 

The  corollary  parameter,  Q^,  the  extinction  efficiency,  defines 
this  relationship  for  the  combined  effects  of  scattering  and  absorptio.1. 

The  extinction  efficiency  is  dependent,  from  Mie  theory,  on  two  parameters, 
the  complex  refractive  index,  and  the  ratio  of  particle  size  to  wavelength 
(r/X)  where  X  is  the  wavelength  of  incident  radiation. 

The  complex  refractive  index  is  commonly  expressed  in  the  form 

m  =  n  +  in ' , 

where  n  and  n*  are  the  real  and  imaginary  parts  respectively.  .  The  real 
part  represents  scattering  while  the  imaginary  part  represents  absorption. 

Volz  (1972)  his  reported  refractive  index  data  for  a  variety  of  aerosol 
substances.  Both  parts  show  considerable  variability  in  their  values 
among  the  different  substances.  In  general,  the  real  part  of  the 
refractive  index  becomes  quite  variable  in  the  middle  infrared  region  and 
increases  slightly  in  the  far  infrared  region.  The  imaginary  part 
increases  markedly  with  increasing  wavelength,  particularly  for  the  water 


WL 


48 


/ 


soluble  continental  aerosols.  Many  aerosol  substances  have  absorption 
indices  in  the  range  0.1  to  0.3  at  8-12  um.  Strongly  absorbing  carbon 
compounds  have  indices  about  0.5.  Even  for  strongly  absorbing  aerosols, 
however,  the  effect  of  refractive  index  variations  on  extinction  efficiency 
is  small  compared  to  the  effect  of  particle  size. 

The  extinction  efficiency  of  aerosol  particles  at  a  specified 
Wavelength  depends  very  strongly  on  particle  size.  This  dependence  is 
usually  expressed  in  terms  of  the  generalized  size  parameter  ,  after  Hanel 
(1976) , 

2irr  .  .  . .  2  .2,1* 

ct£  =  — r—  { (n-1)  +  n*  > 

Fig.  19  shows  the  extinction  efficiency  versus  the  generalized 
size  parameter  for  several  refractive  index  values.  The  value  of  the  size 
parameter  is  primarily  dependent  on  the  ratio  r/A;  we  see  that  the  extinction 
efficiency  is  very  low  for  particles  with  size  parameter  less  than  2. 

Table  8  shows  the  values  of  for  several  particle  sizes  and  refractive 
index  values.  The  table  contains  values  of  extinction  efficiency  taken 
from  Fig.  19.  The  data  show  that  visible  (.55  um)  extinction  efficiency 
is  relatively  uniform  in  the  size  range  indicated.  Extinction  efficiency 
at  10  um, on  the  other  hand,  is  very  low  for  small  particles.  Only  weak 
far  infrared  extinction  occurs  for  particles  smaller  than  2  um,  even  for 
those  with  a  strong  absorption  component  as  represented  by  the  third  value 
of  the  refractive  index  in  Table  8.  This  result  implies  that  when  high 
concentrations  of  the  small  particles  are  present,  even  in  concentrations 
that  may  result  in  rather  law  visibility  values,  the  8-12  u®  transmission 
by  aerosols  may  remain  high. 

The  preceding  discussion  suggests  that  aerosol  extinction  in  the 
8-12  u®  band  depends  primarily  on  the  presence  of  particles  larger  than 
about  2  um.  The  case  studies  presented  in  Section  3  show  that  significant 
numbers  of  particles  in  this  size  range  cam  he  expected  in  the  natural 
aerosol  only  in  the  marine  aerosol  or  when  windblown  dust  particles  are 
present.  These  large  particles  may  also  res  lit  from  man-made  sources  and 
exist  in  the  form  of  dust  or  fly  ash. 
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Fig.  T9.  Efficiency  factor  of  extinction  versus  generalized  size  parameter 
a£  as  computed  from  Mie  theory  for  several  complex  refractive 
indices  n  -  in'.  (Adapted  from  H3nel,  1976). 


The  extinction  efficiency  at  1.06  pm  is  lower  for  particles 
smaller  than  1  um  than  for  the  corresponding  efficiencies  at  0.55  pm. 

The  efficiencies  are  much  higher,  on  the  other  hand,  than  those  at  10  pm 
for  the  same  particles.  Particles  1  pm  or  greater  attenuate  very 
efficiently  at  1.06  pm.  Following  the  argument  above,  these  data  show 
that  attenuation  at  1.06  pm  is  significant  whenever  particles  larger 
than  0.2  pm  are  present  in  significant  concentrations.  This  condition 
is  in  marked  contrast  to  the  requirement  for  particles  larger  than  2  pm 
for  significant  extinction  in  the  8-12  pm  band. 


Table  8 


Size  parameters  (a^)  and  extinction  efficiency  values  (Q£)  as  a 
function  of  refractive  index  (m)  and  particle  radius  at  three  radi¬ 
ation  wavelengths.  The  extinction  efficiencies  were  taken  from 
Fig.  19. 
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6  RELATIVE  HUMIDITY  EFFECTS 


Soluble  aerosol  particles  absorb  liquid  water  on  their  surfaces 
and  grow  in  size  when  relative  humidity  increases.  For  highly  soluble  com¬ 
pounds  when  the  relative  humidity  reaches  a  critical  value,  the  particles 
go  into  liquid  solution  and  a  step  increase  in  particle  size  occurs.  The 
critical  relative  humidity  for  pure  sodium  chloride  is  about  76  percent;  the 
value  for  pure  ammonium  sulphate  is  about  81  percent. 

Particle  growth  with  increasing  relative  humidity  acts  in  two 
ways  to  change  the  aerosol  extinction  properties.  The  principal  effect 
is  through  the  change  in  the  size  distribution.  The  absorption  of  water 
effectively  shifts  the  entire  distribution  of  the  soluble  particles  to 
larger  sizes.  This  increase  in  size  enhances  the  extinction  efficiency 
of  the  particles  at  longer  wavelengths.  The  aggregate  effect  is  signi¬ 
ficantly  greater  extinction  at  those  wavelengths.  A  secondary  but  impor¬ 
tant  effect  of  increasing  relative  humidity  is  the  change  in  the  refrac¬ 
tive  index  of  the  particles.  As  the  particle  water  content  increases, 
the  refractive  index  of  the  particle  approaches  that  of  water. 

This  section  emphasizes  two  important  ideas  regarding  relative 
humidity  effects.  The  first  is  that  the  growth  behavior  of  the  soluble 
particle  depends  on  its  specific  chemical  composition.  The  second  idea 
illustrates  the  effect  of  relative  humidity  increases  in  the  LOWTRAN 
modeled  number  density  and  surface  area  distributions,  compares  these  dis¬ 
tributions  between  the  rural  and  maritime  models,  and  discusses  implica¬ 
tions  for  infrared- extinction. 

■Much  of  the  current  understanding  of  relative  humidity  effects 
derives  from  extensive  theoretical  and  experimental  work  by  Hanel. 

His  early  work,  is  published  in  Volume  19  of  Advances  in  Geophysics  (Hanel, 
1976) .  This  reference  details  the  tueoretical  foundation  for  prediction 
of  particle  growth  with  changing  relative  humidity. 

HSnel  and  Lehmann  (1981)  report  recent  experimental  data  on  the 
dependence  of  growth  on  chemical  composition  jf  the  aerosol.  Figure 
20  presents  a  comparison  of  the  size  dependence  on  relative  humidity  for 
two  different  aerosol  types.  The  curves  were  derived  from  the  observed 
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chemical  composition  of  the  two  aerosols  and  their  predicted  growth  behavior. 
The  Mace  Head,  Ireland  aerosol  was  a  marine  aerosol  with  a  large  sea-salt 
mass  component.  The  growth  curve  shown  is  for  the  sea-salt  particles. 

The  Hohenpeissenberg,  Germany  aerosol  was  a  continental  aerosol  com¬ 
posed  almost  entirely  of  ammonium  sulphate,  ammonium  nitrate  and  sulfuric 
acid. 

The  curves  in  Figure  20  show  that  the  equilibrium  size  of  a 
particle  during  changing  relative  humidity  depends  on  the  direction  of  the 
relative  humidity  change.  This  hysteresis  effect  is  due  primarily  to  the 
existence  of  super-saturated  solutions  in  the  droplet  during  periods  of 
decreasing  relative  humidity.  The  practical  importance  of  the  hysteresis 
effect  to  the  TDA  forecasting  problem  is  small  and  occurs  over  a  narrow 

i 

relative  humidity  range.  The  LOWTRAN  aerosol  models  do  not  account  for 
this  effect 1  Their  tendency  within  this  range  is  to  overestimate  extinc¬ 
tion  slightly  during  increasing  relative  humidity  and  to  underestimate 
during  decreasing  relative  humidity. 

The  primary  feature  of  Fig.  20  is  the  difference  in  growth  behav¬ 
iors  of  the  two  aerosol  compounds.  The  sea-salt  compounds  show  a  much 
stronger  relative  humidity  response  than  the  compounds  in  the  continental 
aerosol. 

The  effect  of  aerosol  response  to  relative  humidity  changes  is 
modeled  in  LOWTRAN  by  changing  the  size  distributions.  Fig.  21  shows 
the  number  density  distributions  of  the  rural  and  maritime  models  at  four 
different  values  of  relative  humidity.  The  total  particle  number  con¬ 
centration  is  fixed  at  15000cm~3  in  the  rural  model  and  at  4000cm  3  in 
the  maritime  model.  In  the  maritime  distribution  one  percent  of  the 
particles  are  of  a  marine  origin. 

Table  9  shows  approximate  values  of  the  number  densities  for  a 
range  of  particle  sizes  for  the  rural  and  maritime  models.  Number  densi¬ 
ties  are  given  at  zero  and  99  percent  relative  humidity  in  each  model. 

The  third  column  under  each  model  shows  the  approximate  ratio  of  number 
density  at  99  percent  to  that  at  zero  percent. 

Note  the  following  features: 

(1)  At  zero  percent  relative  humidity,  the  maritime  model 
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Fig.  20.  Hysteresis  curves  for  maritime  aerosols  (Mace  Head,  Ireland)  and 

continental  aerosols  (Hohenpelssenberg,  Germany) .  Arrows  indicate  the 
direction  of  the  humidity  change  represented  by  each  curve.  (Adapted 
from  Hanel  and  Lehmann,  1981). 

contains  significantly  higher  concentrations  than  the  rural  model  in  the 
size  range  0.8  pm  to  2.0  pm.  The  mode  radius  for  the  sea-salt  particles 
at  zero  percent  relative  humidity  is  approximately  0.2  pm.  The  differ¬ 
ences  in  the  number  concentrations  in  the  range  from  0.8  pra  to  2.0  pm  are 
due  to  the  large  particle  tai-l  of  the  sea-salt  mode  distribution. 

(2)  The  number  densities  in  the  size  range  above  the  2  pm  in¬ 
crease  at  99  percent  relative  humidity  in  much  higher  ratio. in  the  mari¬ 
time  model  than  in  the  rural  model.  This  reflects  the  strong  growth 
response  of  the  sea-salt  aerosol  compared  to  that  of  the  continental 
aerosols.  The  mode  radius  of  the  sea-salt  aerosol  at  99  percent  relative 
humidity  is  about  0.75  pm;  those  of  the  two  modes  in  the  continental 
aerosol  are  about  0.05  pm  and  1.17  pm  respectively.  However,  the  number 
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A  comparison  of  number  density  distributions  in  the  LOWTRAN  rural  and 
maritime  models  at  four  values  of  relative  humidity  (Shettle  and  Fenn 


Table  9 


.  Approximate  number  densities  of  particles  of  various  sizes  at  rela¬ 
tive  humidities  of  zero  percent  and  99  percent  in  the  rural  and  mari¬ 
time  models.  The  third  column  under  each  model  shows  the  approximate 
ratio  (R)  of  the  number  density  at  99  percent  relative  humidity  to 
that  at  zero  percent. 
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concentration  in  the  second  continental  mode  is  much  less  than  the  marine 
aerosol  counts  in  the  maritime  model. 

(3)  The  ratio  of  number  densities  in  the  size  range  below  0.8  pm 
is  somewhat  higher  in  the  rural  model  than  in  the  maritime  model.  This 
difference  reflects  the  dominance  of  the  rural  model  size  distribution  by 
the  continental  aerosol.  The  increase  in  the  number  densities  for  the 
particles  below  6.8  pm  represents  growth  of  the  soluble  continental 
aerosols. 

Referring  back  to  Figures  1  and  2  in  Section  1,  we  may  explain  the 
differences  in  the  extinction  properties  modeled  in  the  rural  and  maritime 
models  as  follows: 
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In  the  8-12  pm  band,  the  strength  of  the  aerosol  extinction 
depends  mostly  on  the  presence  of  the  sea-salt  aerosol.  In  section  5, 
we  showed  that  significant  extinction  Lr.  "his  spectral  region  required 
the  presence  of  particles  larger  than  2  pm.  Table  9  shows  that  these 
particles  are  modeled  in  great  numbers  only  in  the  maritime  model  at  high 
relative  humidities. 

The  relatively  high  numbers  of  the  large  particles  combined  with 
their  high  extinction  efficiency  at  this  wavelength  produce  the  observed 
behavior  of  the  maritime  model's  8-12  pm  transmission  shown  in  Fig.  1. 

At  1.06  pm,  we  observed,  a  significant  difference  in  Fig.  2  in  the 
maritime  model  transmission  compared  to  that  of  the  rural  model  for  a 
given  visibility.  The  relative  humidity  response  of  each  of  the  models 
was  small.  He  showed,  in  Section  5,  that  1.06  pm  extinction  is  dependent 
on  the  number  density  of  particles  larger  than  0.2  pm.  Table  9  shows 
that  these  particles  are  always  present  in  both  models  in  significant 
numbers.  Hence,  the  1.06  pm  extinction  is  much  greater  in  either  model 
than  that  at  8-12  pm. 

The  differences  between  the  maritime  and  rural  model  extinctions 
at  1.06  pm  at  a  given  visibility  are  explained  by  the  presence  or  absence 
of  the  sea-salt  aerosol.  The  size  range  of  the  sea-salt  particles  across 
the  relative  humidity  range  from  zero  to  100  percent  is  such  that  they 
are  effective  attenuators  of  1.06  pm  radiation  at  all  relative  humid¬ 
ities.  Their  presence  in  the  maritime  model  simply  increases  the  re¬ 
lative  concentration  of  effective  1.06  pm  attenuators  at  all  relative 
humidities. 

The  difference  in  the  extinction  properties  of  the  urban  model 
compared  to  the  rural  model  is  of  less  importance  than  the  differences 
due  to  the  marine  aerosol.  The  observed  differeces  between  the  rural 
arid  urban  models  are  due  to  the  soot-like  aerosols  in  the  urban  model. 
Since  these  particles  are  insoluble,  their  effect  is  to  make  the  urban 
model  less  responsive  to  relative  humidity  changes  than  the  rural  model. 
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7  THE  AEROSOL  MODEL  SELECTION  ALGORITHM 


The  Aerosol  Model  Selection  Algorithm  assumes  that  the  most  im¬ 
portant  consideration  in  selection  is  the  presence  or  absence  of  the  sea- 
salt  aerosol.  When  sea  salt  is  thought  to  be  present  in  sufficient  con¬ 
centrations,  then  the  maritime  model  takes  precedence  over  the  rural  and 
urban  models  for  predicting  infrared  extinction,  even  if  urban  pollution 
is  present.  If  sea  salt  is  not  thought  to  be  present  in  sufficient  con¬ 
centrations,  then  the  rural  model  is  re  coemended  for  the  boundary  layer 
except  when  urban  pollutants  are  present;  in  that  case,  the  algorithm 
provides  criteria  for  selection  of  the  urban  model  over  the  rural  model. 

The  criteria  for  determining  whether  sea-salt  aerosols  are  present 
are  based  on  three  considerations:  (1)  the  origin  of  the  air-mass ,  (2)  pos¬ 
sible  addition  of  sea-salt  aerosols  during  aixmass  travel,  and  (3)  eval¬ 
uation  of  mechanisms  for  removal  of  the  sea-salt  aerosol  after  landfall 
of  the  airmass. 

The  role  of  airmass  type  in  defining  the  presence  of  sea  salt  is 
evident  in  data  published  by  Duncan  and  Lindberg  (1981) .  They  fit  ob¬ 
served  particle  size  distributions  in  Europe  to  a  bimodal  distribution 
model  and  performed  Mie  scattering  calculations  at  0.55  pm  and  10  iim. 

Fig.  22,  taken  from  their  report,  presents  a  plot  of  extinction  coeffi¬ 
cients  at  10  pm  versus  extinction  coefficients  at  0.55  pm.  The  size 
distributions  were  measured  in  German  weather  conditions  during  three 
periods:  February-March  1978,  November  1978,  and  February-March  1980. 
Visibility  conditions  varied  from  light  haze  to  dense  fog.  The  data 
represent  measurements  in  three  airmass  types:  continental  polar, 
maritime  arctic,  and  maritime  polar.  Note  the  wide  scatter  of  the  data 
plots  and  the  apparent  separation  by  airmass  type. 

Pigs.  23,  24  and  25  are  plots  of  the  same  data  set  after  the  data 
were  separated  according  to  airmass  type.  The  data  scatter  is  much  re¬ 
duced  in  comparison  with  Pig.  22.  An  overlay  of  these  plots  would  show 
some  differences  between  continental  polar  air  and  the  two  types  of 
maritime  air. 

The  algorithm  assures  addition  of  the  sea-salt  aerosol  when  con¬ 
tinental  airmasses  travel  over  open  ocean  enroute  to  the  point  of  interest 
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Pig.  22. 


Composite  plot  of  extinction  coefficient  at  10  ym  versus  extinction 
coefficient  at  0.55  ym  (units  are  km-1)  (Duncan  and  Lindberg,  1981) . 
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COMPOSITE  PLOT 
MARITIME  ARCTIC  AIR 


EXTINCTION  AT  0.55  yin 


25.  Same  as  Fig.  23  except  for  continental  polar  airaass  (Duncan  and 
Lindberg,  1981). 
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The  length  of  oceanic  trajectory  required  for  use  of  the  maritime  model 
is  assumed  to  be  longer  for  airmasses  traveling  over  cold  water  than  for 
those  traveling  over  warm  water.  The  specific  values  of  800  km  and  500  km 
for  trajectory  length  are  only  estimates.  This  criterion  is  a  means  of 
accounting  for  the  effect  of  static  stability  on  the  relative  effective¬ 
ness  of  mechanisms  for  vertical  transport  of  the  sea-salt  particles. 

Two  mechanisms  for  removal  are  considered.  The  algorithm  assumes 
that  most  sea-salt  particles  larger  then  2  iim  will  be  lost  by  sedimenta¬ 
tion  within  three  days  after  landfall  of  the  airmass.  The  algorithm  also 
assumes  that  any  precipitation  of  moderate  or  greater  intensity  will  re¬ 
move  most  of  the  sea-salt  particles.  If  either  of  the  above  criteria  is 
met,  then  the  rural  model  is  recommended,  even  for  airmasses  of  maritime 
origin. 

Urban  model  selection  criteria  are  based  on  the  strength  of  the 
urban  pollution  source,  the  effectiveness  of  ventilation  by  the  boundary 
layer  wind,  and  effect  of  static  stability  on  vertical  dispersion.  Ap¬ 
pendix  A  describes  how  the  selection  criteria  were  developed. 

Fig.  26  presents  the  logic  chart  which  constitutes  the  selection 
algorithm.  The  following  briefly  describes  the  organization  of  the  chart. 

(a)  Fig.  26-A  separates  airmasses  by  origin,  and  treats  the 
possible  transformation  of  an  airmass  with  a  continental  origin  and  an 
oceanic  trajectory. 

(b)  Fig.  26-B  treats  the  mechanisms  for  removal  of  the  sea-salt 
aerosol  from  a  maritime  airmass,  namely,  sedimentation  and  washout. 

(c)  Figs.  26-C  and  26-0  treat  the  problem  of  determining  when 
the  urban  model  should  be  used  to  describe  a  polluted  continental 
aerosol.  The  urban  model  logic  arrives  at  a  critical  value  of  the  fore¬ 
cast  boundary  wind  speed.  When  the  forecast  wind  speed  is  less  than  the 
critical  valuer  use  of  the  urban  model. is  recommended. 


Fig.  26.  The  aerosol  model  selection  process 
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Fig.  26.  (Continued) 
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8  SUMMARY  AND  CONCLUSION 

Section  1  describes  how  the  LOW-’.’RAN  aerosol  models  show  the  depend¬ 
ence  of  infrared  extinction  on  the  nature  of  the  aerosol.  Effective  ap¬ 
plication  of  LOWTRAN  in  operational  prediction  of  infrared  extinction 
requires  some  means  of  estimating  aerosol  properties  to  aid  in  "best” 
aerosol  model  selection. 

In  Section  2  we  described  various  methods  of  characterizing  aer¬ 
osols.  The  total  number  concentration  and  number  density  versus  size 
distribution  has  a  primary  effect  on  the  spectral  dependence  of  extinction. 
Hater  solubility  and  specific  chemical  composition  are  the  primary  factors 
in  determining  the'  aerosol  growth  response  to  changes  in  relative  humidity. 

The  descriptive  case  studies  in  Section  3  show  the  principle  dif¬ 
ferences  in  various  types  of  aerosols.  The  sudies  showed  the  existence 
of  the  global  background  aerosol.  This  component  of  the  aerosol  shows 
little  variability  in  space  and  time.  The  variable  aerosol  components 
are  superimposed  on  the  background  aerosol  and  consist  principally  of 
three  aerosols:  sea  salt,  urban  pollutants  and  windblown  dust. 

Section  4  describes  the  processes  that  determine  the  important 
time  and  space  variability  in  aerosol  properties.  Continuous  steady-state 
production  mechanisms  maintain  the  global  background  aerosol.  Sea-salt 
aerosol  production  by  the  wind,  aerosol  production  in  urban  and  industrial 
complexes,  and  windblown  dust  are  the  most  important  sources  of  the  variable 
components .  The  principle  aerosol  removal  mechanisms  on  a  large  scale 
are  sedimentation  and  rainout/washout.  .  Dispersion  is  important  in  dilution 
of  urban  aerosol  concentrations.  We  found  that  precipitation  and  washout 
are  most  effective  in  removal  of  the  aerosols  larger  than  2  ym. 

In  Section  5  we  showed  the  importance  of  the  number  density  of 
various  particle  sizes  in  determining  the  wavelength  dependence  of  extinc¬ 
tion.  Strong  aerosol  extinction  at  8-12  ym  requires  significant  concen¬ 
trations  of  particles  larger  them  2  ym.  Near  IR  extinction  (e.g.,  1.06 
ym)  depends  primarily  on  the  concentration  of  particles  in  the  size  range 
from  0.2  ym  to  0.8  ym. 

Section  6  emphasized  two  important  points  regarding  aerosol 
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growth  with  increasing  relative  humidity.  The  growth  behavior  of  soluble 
aerosols  depends  on  their  chemical  composition;  sea-salt  aerosols  exhibit 
a  growth  response  much  stronger  than  that  of  continental  aerosols.  The 
change  in  number  density  distribution  with  increasing  relative  humidity 
results  in  a  strong  relative  humidity  dependence  of  8-12  pm  extinction  in 
marine  aerosols. 

Relative  humidity  dependence  at  1.06  pm  on  the  other  hand  is  rela¬ 
tively  weak.  The  maritime  aerosols  exhibit  much  stronger  extinction  than 
the  continental  aerosols.  This  difference  is  explained  by  the  presence 
of  sea  salt,  which  is  an  effective  attenuator  at  1.06  pm  across  its  entire 
range  of  particle  sizes. 

The  general  conclusion  is  that  the  primary  contributor  to  varia¬ 
bility  in  infrared  extinction  is  the  sea-salt  aerosol .  The  first-order 
forecasting  problem  is  to  assess  the  presence  of  this  component.  Studies 
of  known  production  and  los3  mechanisms  and  of  the  dependence  of  extinc¬ 
tion  properties  pn  airmass  type  indicate  that  a  reasonable  assessment  can 
be  made  from  an  analysis  of  airmass  history.  Airmass  analysis  is  the 
primary  basis  of  the  Aerosol  Model  Selection  Algorithm. 
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APPENDIX  A  DEVELOPMENT  OF  CRITERIA  FOR  SELECTION  OF  THE  URBAN  AEROSOL  MODEL 


Diffusion  theory  states  that  the  concentration  of  aerosols  in  and' 
near  an  urban  area  is  a  complex  function  of  the  strength  of  the  pollution 
sources,  the  horizontal  diffusive  properties  of  the  atmosphere,  the  wind 
speed  and  the  vertical  static  stability.  Pasquill  (1961,  1971)  developed 
the  early  theory  and  simple  graphical  techniques  for  evaluation  of  down¬ 
wind  concentrations  from  a  point  source.  The  well-known  Pasquill  stability 
classifications  are  a  part  of  his  solution  to  the  point  source  problem. 
Other  authors  (McElroy,  1969,  Turner,  1964)  amplified  on  his  original 
development  and  developed  more  complex  models  of  pollutant  concentrations 
from  an  array  of  point  sources.  These  models  predict  urban  aerosol 
concentrations  id  large  urban  complexes. 

Studies  of  measured  pollution  concentrations  (Zib,  1980;  Gifford 
and  Hanna,  1973;  Hanna,  1978)  in  large  urban  complexes  have  shown  that 
the  urban  aerosol  concentration  may  be  estimated  with  reasonable  accuracy 
using  the  simple  formula: 

CO 

X  *  V  (A-l) 

where  x  is  the  concentration  value  within  the  urban  complex,  C  is  a 
dimensionless  parameter ,  Q  is  the  pollutant  source  strength  in  units  of 
mass  per  unit  area  per  unit  time,  and  u  is  the  boundary  layer  wind  speed. 

An  average  value  of  C  *  225  has  been  found  to  work  well  overall.  The 
parameter  C  has  been  found  to  depend  on  the  size  of  the  urban  area  and 
on  static  stability  in  the  boundary  layer. 

Gifford  and  Hanna  (1973)  defined  an  expression  for  C  which  depends 
on  the  distance  of  a  point  from  the  upwind  edge  of  the  pollutant  source 
area.  Hanna  (1978)  used  this  definition  in  a  theoretical  analysis  of 
this  parameter  and  its  dependence  on  static  stability.  The  parameter  D 
is  the  distance  from  a  point  within  the  urban  complex  to  its  upwind  edge 
In  a  large  city,  the  expected  value  is  one-half  the  characteristic  diameter 
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Table  A-l .  Calculated  values  of  C  for  a  range  of  stability  parameters  and 

city  sizes,  D.  The  parameters  a  and  b  are  used  to  model  vertical 
dispersion.  (Adapted  from  Hanna,  1978.)  . 


Meteorological 

conditions 

a 

b 

C  for 

D  *  10  m 

C  for 

D  -  10  m 

C  for 

D  =  105m 

Very  unstable 

0.40 

0.91 

41 

51 

63 

Unstable 

0.33 

0.86 

46 

63 

87 

Neutral 

0.22 

0.80 

73 

115 

182 

Stable 

0.06 

0.71 

341 

662 

1301 

of  the  source  area.  He  used  the  relation 

a  -  axb  (A- 2) 

z 

where  a  ,  the  vertical  dispersion  parameter ,  represents  thi  effects  of 
z 

static  stability  on  C.  The  x  coordinate  is  measured  in  the  direction  of 
wind  from  the  upwind  edge  of  the  source  area.  Table  A-l  presents  the 
results. 


Hanna  used  reported  values  of  suspended  particle  data  for  several 
cities  to  verify  this  theoretical  values.  He  found  an  overall  average 
value  of  C  *  200  to  be  in  good  agreement.  He  also  derived  diurnal  curves 
of  C  from  observed  values  of  carbon  monoxide  concentration  in  New  Jersey, 
Maryland  and  Colorado.  He  showed  that  mpst  of  the  diurnal  variation  in  C 
is  due  to  variation,  in  static  stability.  He  used  corresponding  curves  of 
the  diurnal  variation  in  wind  speed  averaged  for  Newark,  Philadelphia, 

Denver  to  calculate  the  C  curves.  He  tested  these  curves  by 
calculating  concentrations  of  carbon  monoxide  in  Los  Angeles  and  comparing 
them  to  measurements.  His  results  were  as  accurate  overall  as  those  from 
a  much  more  ©  implex  computer  simulation  model. 

Gifford  and  Hanna  (1973)  published  pollutant  concentration  data 
for  a  number  of  American  urban  areas  of  various  sizes  and  degrees  of  in¬ 
dustrialization.  We  extracted  a  subset  of  these  data  in  search  of  a  simple 
method  by  whioh  the  user  of  the  LOWTRAN  aerosol  models  may  evaluate  the 
of  urban  aerosols  in  an  urban  complex  at  a  specified  time. 


concentration 


Table  A-2  presents  the  data  used  in  this  analysis. 
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He  computed  the  values  of  the  source  strength  per  unit  area  using 
the  total  source  strength  for  particles  and  SO^  combined  and  the  areas  pro¬ 
vided  by  the  Department  of  Health,  Education  and  Welfare  (DHEW)  as  shown 
in  Column  1.  The  resulting  values  of  source  strength  per  unit  area  are 
shown  in  Column  5 . 

Pig.  A-l  is  a  plot  of  the  average  annual  concentrations  of  particles 
and  SO^  against  the  average  source  strength  per  unit  area  for  each  of  the 
12  urban  areas  listed  in  Table  A-2-  City  names  and  area  size  of  the  pol¬ 
lutant  sources  (Column  1)  are  indicated  at  each  data  point.  Data  for  all 
but  four  of  the  urban  areas  approximate  a  linear  relationship..  The  four 
outliers  are  values  for  metropolitan  areas  with  heavy  industrialization 
and  very  large  areas  of  pollutant  sources.  The  two  dashed  lines  in  Pig. 

A-l  are  drawn  to  approximate  the  relationship  between  pollutant  concentra¬ 
tions  and  source  strength.  The  lower  line  represents  the  relationship  for 
small  to  medium  sized  urban  areas.  The  upper  dashed  line  represents  the 
relationship  for  large  heavily  industrialized  urban  complexes. 

These  dashed  lines  are  represented  by  the  linear  equations  of  the 

form: 


X  -  a  +  b Q 


(A-3) 


X  m  230  +  6.2  Q  for  large  urban  areas  (A-4) 

X  “  118  4.6  Q  for  small  to  medium  size  cities.  (A-5) 

The  slope  corresponds  to  the  value  of  C/U  in  Eq.  A-l.  Using  the 
average  annual  wind  speed  of  7.6  m  sec  ^  for  large  cities  and  7.1  n  sec  * 
for  small-medium  sized  cities,  the  resulting  values  of  C  are  47  and  32 
respectively. 

Eq.  (A-l)  requires  that  the  linear  relationship  between  concen¬ 
tration  and  source  strength  pass  through  the  origin.  Eq.  (A-4)  and  (A-5), 
on  the  other  hand,  have  a  large  y  intercept.  The  resulting  value  of  C  is 
much  less  than  that'  found  when  the  form  in  Eq.  (A-l)  is  used. 

We  assumed  that  values  of  C  from  Eqs.  (A-4)  and  (A-5) ,  were  re¬ 
presentative  of  neutral  stability  conditions.  Furthermore,  we  assumed 
that  an  average  C  value  of  40  would  be  representative  for  cities  of  all 
sizes  under  neutral  stability.  Using  the  data  by  Hanna  in  Table  A-l, 
we  estimated  values  of  C  ■  20  for  unstable  conditions  and  C  *  120  for  stable 
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•  Minn-St.  Paul  (775) 


Table  A-2.  Source  data  for  12  cities  used  in  evaluating  the  role  of  source 

strength  and  size  of  the  urban  area.  Data  were  taken  from  Gifford 
and  Hanna  (1973) . 


Providence  155  141  238  8.3  26.1 
Denver  260  58  135  6.3  6.4 
Chicago  2590  2366  360  7.3  2 u_3 
St.  Louis  595  ,  838  293  6.5  40.5 
Philadelphia  4400  1399  341  7.8  9.1 


(1) 


2 

Approximate  area  (km  )  enclosed  by  an  urban  particulate  source  with 
strength  greater  than  0.1  tons  day-'*'  mi-^. 


(2) 


Total  average  annual  strength  of  the  urban  pollution  sources  (103 4  tons 
yr_i) .  This  includes  both  particles  and  SO^  emissions. 


(3)  Total  average  annual  concentrations 

tides  and  SO  concentrations, 
x 

(4)  Average,  annual  wind  speed  (m  sec*'3'). 


Q  for  particles  and  SO  (yg  m" 


This  value  includes  par- 


conditions  as  representative  values  for  use  in  Eq.  (A-3)  for  cities  of  all 
sizes.  Since  b  =  C/U,  we  may  rewrite  this  equation  in  the  form 


If  we  assume  that  criteria  for  selection  of  the  urban  aerosol 
model  in  lieu  of  the  rural  model  may  be  based  on  some  critical  value  of 
the  concentration,  then  the  selection  criteria  may  also  be  expressed  in 
terms  of  the  boundary  layer  wind  speed  if  the  stability  category  is  de¬ 
fined. 

We  found  little  data  on  the  values  of  pollutant  concentration  at 
which  the  urban  aerosol  model  should  be  used.  Junge  (1963)  cited  typical 
values  of  200-800  yg  m  ^  under  moderate  pollution  conditions.  He  reported 
concentrations  as  high  as  4000  yg  m  ^  in  severe  London  smog.  These  compare 
to  typical  concentrations  of  5-50  yg  m  ^  for  the  rural  background  aerosol. 
We  selected  a  critical  value  of  400  yg  m  ^  as  criteria  for  selection  of 
the  urban  aerosol  model.  Using  Eq.  (A-6) ,  the  values  of  C  described  above, 
and. the  values  of  a  in  Eqs.  (A-4)  and  (A-5) ,  we  solved  for  a  critical 
value  of  boundary  layer  wind  speed  for  each  stability  condition  for  large 
cities  and  for  small  to  medium  size  cities.  The  resultant  wind  speeds 
provide  guidance  for  selecting  the  urban  aerosol  model  on  the  basis  of 
three  parameters,  namely,  the  size  and  degree  of  industrialization  of 
the  urban  area,  the  stability  category  and  the  forecast  boundary  layer 
wind  speed.  The  critical  wind  speeds  are  presented  in  Fig.  26-0  in  the 
main  body  of  the  report. 


